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Abstract 

We study the diffusion front for a natural two-dimensional model where 
many particles starting at the origin diffuse independently. It turns out 
that this model can be described using properties of near-critical perco- 
lation, and provides a natural example where critical fractal geometries 
spontaneously arise. 

1 Introduction 

In this paper, we study a simple two-dimensional model on a planar lattice where 
a large number of particles, all starting at a given site, diffuse independently. 
As the particles evolve, a concentration gradient appears, and the random inter- 
faces that arise can be described by comparing the model to an inhomogeneous 
percolation model (where the probability that a site is occupied or vacant de- 
pends on its location). We exhibit in particular a regime where the (properly 
rescaled) interfaces are fractal with dimension 7/4: this model thus provides 
a natural setting where a (stochastically-) fractal geometry spontaneously ap- 
pears, as observed by statistical physicists, both on numerical simulations and 
in real-life situations. 

To our knowledge, the study of the geometry of such "diffusion fronts" has 
been initiated by the physicists Gouyet, Rosso and Sapoval in [20]. In that 
1985 paper, they showed numerical evidence that such interfaces are fractal, 
and they measured the dimension Df = 1.76 ± 0.02. To carry on simulations, 
they used the approximation that the status of the different sites (occupied 
or not) are independent of each other: they computed the probability p(z) of 
presence of a particle at site z, and introduced the gradient percolation model, an 
inhomogeneous percolation process with occupation parameter p(z). Gradient 
percolation provides a model for phenomena like diffusion or chemical etching 
- the interface created by welding two pieces of metal for instance - and it has 
since then been applied to a wide variety of situations (for more details, see the 
recent account [7] and the references therein). 

Based on the observations that the front remains localized where p{z) ~ p c 
and that the dimension Df was close to 7/4, a value already observed for critical 
(homogeneous) percolation interfaces, they argued that there should be a close 



similarity between these diffusion fronts and the boundary of a large percolation 
cluster. For critical percolation, this conjectured dimension 7/4 was justified 
later via theoretical physics arguments (see for instance [HO [2]). In 2001, 
thanks to Smirnov's proof of conformal invariance [H] and to SLE computations 
by Lawler, Schramm and Werner (see e.g. [321 [H]), a complete rigorous proof 
appeared [3] in the case of the triangular lattice. It is actually known that in 
this case, the interfaces are described by the so-called SLE(6) process, a process 
which has Hausdorff dimension 7/4. There is also a discrete counterpart of this 
result: the length of a discrete interface follows a power law with exponent 7/4, 
that comes from the arm exponents derived in [22]. For more details, the reader 
can consult the references \W\ [24] . 

In the paper [15J, we have studied the gradient percolation model from a 
mathematical perspective, building on the works of Lawler, Schramm, Smirnov 
and Werner. Combining their results with Kesten's scaling relations [9], one first 
gets a rather precise description of homogeneous percolation near criticality in 
two dimensions |22II16|. for instance the power law for the so-called characteristic 
length 

£b)-b-i/2r 4/3+o(l) 

as p — ► 1/2 (here and in the sequel we restrict ourselves to the triangular lattice 
T). Estimates of this type (and the underlying techniques) then enabled us to 
verify the description in [20J for gradient percolation, and to confirm rigorously 
that dimension 7/4 shows up in this context. 

The goal of the present paper is to investigate "diffusion front" models math- 
ematically, and to use our results concerning gradient percolation to show that 
their geometry and roughness can also be rigorously described. 

Let us now present our main model in more detail. We start at time t = 
with a large number n of particles located at the origin, and we let them 
perform independent (simple) random walks. At each time t, we then look at 
the sites containing at least one particle. These occupied sites can be grouped 
into connected components, or "clusters", by connecting two occupied sites if 
they are adjacent on the lattice. 

As time t increases, different regimes arise, as depicted in Figure [TJ Roughly 
speaking and in the large n-limit, at first, a very dense cluster around the 
origin forms. This clusters grows as long as t remains small compared to n. 
When t gets comparable to n, the cluster first continues to grow up to some 
time t max = A max n, then it starts to decrease, and it finally dislocates at some 
critical time t c — X c n - and it never re-appears. The constants A max < A c can 
be explicitly computed, and even if they depend on the underlying lattice, their 
ratio is (expected to be) universal, equal to l/e. 

To be more specific, the phase transition corresponding to A c can be de- 
scribed as follows: 

• When t = An for A < A c , then (with probability going to one as n — > oo), 
the origin belongs to a macroscopic cluster. Its boundary is roughly a 
circle - we will provide a formula for its radius in terms of A and n, it is 
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(a) t = 10. 



(b) t = 100. 



(c) t = 500. 




(d) t = 1000. 



(e) t = 1463 (= |AmaxnJ). 



(f) t = 2500. 




(g) t = 3977 (= LAcnJ). 



(h) t = 5000. 



(i) t = 10000. 



Figure 1: Different stages of evolution as the time t increases, for n = 10000 
particles. 
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of order ^fn - and its roughness can be described via the dimension 7/4. 
This is a dense phase. 

• When t = An for A > A c , then all the clusters - in particular the cluster 
of the origin - are very small, and the whole picture is stochastically 
dominated by a subcritical percolation picture. This is a dilute phase. 

Actually, alongside these two principal regimes, at least two other interesting 
cases could be distinguished: the (near-)critical regime when t is very close to 
X c n, and a short-time transitory regime when t is very small compared to n 
(t <C n a for any a > 0, typically t = logn). We do not discuss these regimes in 
the present paper. 

We finally study a variant of this model, where now a source creates new 
particles at the origin in a stationary way - one can for instance think of the 
formation of an ink stain. If the new particles arrive at a fixed rate fi > 0, then 
after some time, one observes a macroscopic cluster near the origin, that grows 
as time passes. We explain briefly how such a model can be described using the 
same techniques and ideas. 

To sum up, our results are thus strong evidence of universality for random 
shapes in the "real world". They indeed indicate that as soon as a density gra- 
dient arises, random shapes similar to SLE(6) should appear - at least under 
some hypotheses of regularity for the density, and of (approximate) spatial in- 
dependence. However, note that it is explained in [T7l [18] that these random 
shapes are also somewhat different from SLE(6): they are locally asymmetric. 

2 Setup and statement of the results 
2.1 Description of the model 

As previously mentioned, our results will be obtained by comparing our model to 
an (inhomogeneous) independent percolation process, for some well-chosen oc- 
cupation parameters. This process is now rather well-understood, which allows 
one to get precise estimates on macroscopic quantities related to the interface, 
such as the size of its fluctuations or its length. In particular, the interface is 
a fractal object of dimension 7/4, whose geometry is close to the geometry of 
SLE(6), the Schramm-Loewner Evolution with parameter 6. 

Recall that site percolation on a lattice G = {V,E), where V is the set of 
vertices (or "sites"), and E is the set of edges (or "bonds") can be described as 
follows. We consider a family of parameters (p(z), z € V), and we declare each 
site z to be occupied with probability p(z), vacant with probability 1 — p{z), 
independently of the other sites. We are then often interested in the connectivity 
properties of the set of occupied sites, that we group into connected components, 
or "clusters". In the following, we use the term "interface" for a curve on the 
dual hexagonal lattice bordered by occupied sites on one side, and by vacant 
sites on the other side (the boundary of a finite cluster for instance). 
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In the standard homogeneous case where p(z) = p 6 [0, 1], it is known for 
a large class of lattices that there is a phase transition at a certain critical 
parameter p c e (0, 1) (depending on the lattice), and it is a celebrated result of 
Kesten [S] that p c = 1/2 on the triangular lattice T = (V T ,E T ): when p < 1/2 
there is (a.s.) no infinite cluster, while for p > 1/2 there is (a.s.) a unique 
infinite cluster. At the transition exactly, for the value p = p c of the parameter, 
lies the critical regime. This regime turns out to be conformally invariant in 
the scaling limit, a very strong property established by Smirnov [21] that allows 
one to prove that cluster boundaries converge to SLE(6). One can then use 
SLE computations O Q3] to derive most critical exponents associated with 
percolation |22j . 

From now on, we focus on the lattice T, since at present this is the lattice for 
which the most precise properties, like the value of the critical exponents, are 
known rigorously - a key step being the previously-mentioned conformal invari- 
ance property [21], a proof of which is still missing on other lattices. However, 
let us mention that most results remain valid on other common regular lattices 
like Z 2 , except those explicitly referring to the value of critical exponents, like 
the fractal dimension of the boundary (see Section 8.1 in [16| for a discussion of 
this issue). 

As usual in the statistical physics literature, the notation / xj means that 
there exist two constants C\,Ci > such that C\g < f < C25, while / « g 
means that log //logy — ► 1 - we also sometimes use the notation which 
has no precise mathematical meaning. 

Let us now describe the model itself. We assume that at time t = 0, a fixed 
number n of particles are deposited at the origin, and we let these particles 
perform independent random walks on the lattice, with equal transition proba- 
bilities 1/6 on every edge. We then look, at each time t, at the "occupied" sites 
containing at least one particle - note that we do not make any exclusion hy- 
pothesis, a given site can contain several particles at the same time. We denote 
by ir t = (vrt(^), z e V T ) the distribution after t steps of a simple random walk 
starting from 0, so that the probability of occupation for a site z is 

p n , t {z) = 1 - (1 - 7T t (z)) n ~ 1 - e- n ^ z \ (2.1) 

and we expect 

/ \ Ct -\\ z \\ 2 /t 

for some constant Ct depending only on the lattice. For our purpose, we will 
need a rather strong version of this Local Central Limit Theorem type estimate, 
that we state in the next section. We also denote by N n> t(z) the number of 
particles at site z, at time t. 

Roughly speaking, p n ,t(z) is maximal for z = 0, and p n ,t(0) decreases to as 
the particles evolve. Hence, different regimes arise according to t: if p n ,t(0) > 
1/2, the picture near the origin is similar to super-critical percolation. There 
is thus a giant component around the origin, whose boundary is located in a 
region where the probability of occupation is close to p c — 1/2. On the other 
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hand, if p n ,t(0) < 1/2 the model behaves like sub-critical percolation, so that 
all the different clusters are very small. 

2.2 Main ingredients 

Local Central Limit Theorem 

To describe the fine local structure of the boundary, we need precise estimates 
for 7r t . We use the following form of the Local Central Limit Theorem: 



(2.2) 



uniformly for z G V T , ||z| < < 9 / 16 . This equation means that locally, the 
estimate given by the Central Limit Theorem holds uniformly: the constant 
V3/2 comes from the "density" of sites on the triangular lattice. 

This result is sharper than the estimates stated in the classic references 



(231 HI], in particular for atypical values like Cy/t logf, and such a strong form 
is needed below. It can be found in Chapter 2 of [12j, where estimates are derived 
for random walks with increment distribution having an exponential moment. 
In that book, results are stated for Z 2 , but they can easily be translated into 
results for the triangular lattice T by considering a linear transformation of T 
onto Z 2 . One can take for instance the mapping x + e l7r ^y i— > x + iy (in complex 
notation), which amounts to embedding T onto Z 2 by adding NW-SE diagonals 
on every face of Z 2 - one then has to consider the random walk with steps 
uniform on {(±1, 0), (0, ±1), (±1, qpl)}, 

We also use several times the following a-priori upper bound on 7r t (z), pro- 
viding exponential decay for sites z at distance > i 1 / 2 ^: for some universal 
constant C > 0, we have 

7rt(z) < Ce~^ 2 ' 2t (2.3) 

for all sites z, for any time t > 0. This bound is for instance a direct consequence 
of Hoeffding's inequality. 



Poisson approximation 

Note that if we take a Poissonian number of particles N ~ V(n), instead of 
a fixed number, the thinning property of the Poisson distribution immediately 
implies that the random variables (N ni t(z),z G V T ) are independent of each 
other, with N nit (z) ~ V(rnr t (z)). We thus simply get an independent percola- 
tion process with parameters 

Pn ,t{z) = l-e- n "*W. (2.4) 

As will soon become clear, this Poisson approximation turns out to be valid, 
since we only need to observe a negligible fraction of the particles to fully de- 
scribe the boundary: what we observe locally looks very similar to the 'Poisso- 
nian case". For that, we use the classic Chen-Stein bound (see Theorem 2.M on 
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p. 34 in [4], or [6]). Assume that W = X\ + . . . + X^ is a sum of fc independent 
random variables, with Xi ~ B(pi). Then for A = E[VT] = X^=iP»> we have for 
the total variation norm 

\\C(W) - P(A)|| ry < l —^j^ P l (2.5) 

i— 1 

Now, assume that we look at what happens in a (small) subset A of the sites. 
The distribution in A can be obtained by first drawing how many particles fall 
in A, and then choosing accordingly the configuration in A. Each particle has a 
probability n t (A) = J^zeA^ti 2 ) to ^ e m A so tnat the number of particles in 
A is iVt(.A) ~ B(n,ir t (A)), and the Chen-Stein bound implies that (A = nwt(A) 
in this case) 

\\C(N t (A)) - V(nn t (A))\\ TV < —L^mr t (A) 2 = n t (A). (2.6) 

We can thus couple N t (A) with ~ V(mr t (A)) so that they are equal with 
probability at least 1 — 7r t (A), which tends to 1 if A is indeed negligible. Hence 
to summarize, we can ensure that the configuration in A coincides with the 
"Poissonian configuration", i.e. with the result of an independent percolation 
process, with probability at least 1 — ir t (A). 



Radial gradient percolation 

We use earlier results on near-critical percolation, especially Kesten's paper [9]. 
All the results that we will use here are stated and derived in [16], and we 
follow the notations of that paper. In particular, we use the basis (l,e i7r / 3 ), 
so that [01,02] x [61,62] refers to the parallelogram with corners a, + bje l7T ^ 3 
{1,2}). 

We recall that the characteristic length is defined as 

L e (p) = min{n s.t. F p (C H ([0,n] x [0,n])) < e} (2.7) 

when p < 1/2 (Ch denoting the existence of a horizontal occupied crossing), 
and the same with vacant crossings when p > 1/2 (so that L t {p) = L e (l — p)), 
for any e € (0, 1/2), and that we have 

L e (p)~LAp)*\p-l/2r 4/3 (2-8) 

for any two fixed < e,e' < 1/2. We will also use the property of exponential 
decay with respect to L: for any e g (0, 1/2) and any k > 1, there exist some 
universal constants d = d(k,e) such that 

P P (Cjj([0, n] x [0, kn])) < C x e' c ^ n ^^ (2.9) 

uniformly for p < 1/2 and n > 1, In the following, we fix a value of e, for 
instance e = 1/4, and we simply write L, forgetting about the dependence on e. 
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The aforementioned approximation allows to use for our model results con- 
cerning the independent Bernoulli case: we hence consider a "radial" gradi- 
ent percolation, a percolation process where the parameter decreases with the 
distance to the origin. Let us first recall the description obtained in a strip 
[0,^at] x [0,N] with a constant vertical gradient, i.e. with parameter p(y) = 
1 — y/N: we will explain later how to adapt these results in our setting. The 
characteristic phenomenon is the appearance of a unique "front", an interface 
touching at the same time the occupied cluster at the bottom of the strip and 
the vacant cluster at the top: if In > N 4 / 7+s for some 5 > 0, this happens with 
high probability (tending to 1 as n — > oo). Moreover, for any fixed small e > 0, 
also with high probability, 

(i) this front remains localized around the line y = y c = \_N/2\ , where p(y) ~ 
1/2, and its fluctuations compared to this line are of order 

a N = sup {a s.t. L(p(y c + a)) > a) « N 4/1 \ (2.10) 

more precisely they are smaller than N 4 ^ 7+€ (the front remains in a strip 
of width 2iV 4 / 7+£ around y — y c ), and larger than 7V 4 / 7-£ (both above 
and below y — y c ), 

(ii) and since it remains around the line y = y C) its behavior is comparable 
to that of critical percolation, so that in particular its fine structure can 
be described via the exponents for critical percolation: for instance, its 
discrete length Lm satisfies 

N 3/7 - e £ N <L N < N 3/7+£ £ N . (2.11) 



2.3 Statement of the results 

In view of the previous sections, let us introduce 

*.(r) = £ t e- r2/t (2.12) 

(so that 7r t (z) = 7f t (||z||)[l + (9(l/i 3 / 4 )] uniformly for ||z|| < t 9 ' 16 by Eq.|H|), 
and the associated "Poissonian" parameters 

p n ,t(r) = l-e- n ^ r l (2.13) 

Since p n ^ is a decreasing continuous function of r, and tends to when 
r — > oo, there is a unique r such that p n ,t(r) = 1/2 iff p n< t(0) > 1/2, or 
equivalently iff t <t c — X c n, with 

X C = ^T 2 (- 0.397...). (2.14) 

We introduce in this case 

< t = (p n ,t)- 1 (l/2), (2.15) 



S 



where (p n ,t) 1 denotes the inverse function of p n ,t- In particular, r* t = when 
t = X c n. One can easily check that 



<>< = f l0S ^' (2 " 16) 

and a small calculation shows that r* t , as a function oft, increases on (0, A max n] 
and then decreases on [A max n, A c n] (finally being equal to 0), with 

A max = ^ = -^— (-0.146...)- (2.17) 
e 2e7rlog2 

In particular, the ratio i max /t c is equal to 1/e, and is thus independent of n - 
and though the lattice appears in A c and A max , as the "density" of sites, V3/2 
here, it disappears when we take the ratio t max /t c . For future reference, note 
also that 

• when t — An, for some fixed A < A c , then r* ( increases as \fi x sjri: 



log y \ft, 
when t = n a , for some < a < 1, 

Vt 



ViJ - ij logt + logA c = - lv/tlogt + O^ 



We also introduce 

<T^ t = sup {a s.t. L(p(Lr;, t J ± a)) > o) (2.18) 

that will measure the fluctuations of the interface. We will see that a^ t w i 2 / 7 
under suitable hypotheses on n and t. We are now in a position to state our 
main result. We denote by S r (resp. S r , r ') the circle of radius r (resp. the 
annulus of radii r < r') centered at the origin. 

Theorem 1. Consider some sequence t n —> oo. Then 

(i) If t n > An /or some A > A c , £/ien there exists a constant c = c(A) smc/i 
£/ia£ 

Wn,t{every cluster is of diameter < clogn) — ► 1 

as n — > oo. 

fii) If e ( l °s n ) a <t n < A'n /or some a G (0, 1) and A' < A c , £/ien /or any /ixed 
smaZZ e > 0, ifte following properties hold with probability tending to 1 as 

n — > oo; 

fa,) £/iere is a unique interface j n surrounding dSt* , 
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(b ) this interface remains localized in the annulus of width 2t n ' £ around 
9S r * (recall that r* { ~ \ft^), and compared to this circle, its 

2 II— € 

fluctuations, both inward and outward, are larger than t n , 

(c) its length - number of edges - L n satisfies t^/ 7 e < L n < tf/ 7+€ (in 
other words, it behaves as L n rj t 5 J~' ' ). 



3 Proofs 

From now on, we choose not to mention the dependence of t on n in most 
instances, for notational convenience. 



Case t n > Xn for some A > A c 

Let us start with case ©, which is the simpler one. In this case, the hypothesis 
t n > Xn implies that for any site z, 

Pn,t(\\4)<Pn,t(0) < l-e-^ log2 =p< 1/2. (3.1) 
For z such that \z\ < t 9 ^ 16 , since 

n t (z) = n(\\4)[i + 0(iA 3 / 4 )] - n t (\\4) + o(i/« 7/4 ) 

uniformly, we also have 

Pn j(z) = 1 - e - n * t{z) = 1 - e -«ff t (ll*ll)+o(i/« 3/4 ) < p > < 1/2 (3.2) 

for some p' > p. Moreover, this bound also holds for ||z|| > t 9 / 16 by Eg. (|2.3|) . 

This enables us to dominate our model by a subcritical percolation process. 
Take some S > 0, and consider the model with N ~ V((l + S)n) particles: 
N > n with probability tending to 1, and in this case the resulting configuration 
dominates the original one. Moreover, if we choose 5 such that 

A c (l + S)n<^±±n, 

the hypothesis t n > A(l + S)n holds for some A € (A c , A), so that the bound 
l|3.2p is satisfied by the associated parameters for some p" G (p' ,1/2). 

Let us assume furthermore that t < n 3 , and subdivide St into ~ t 2 horizontal 
and vertical parallelograms of size f logn x clogn. Any cluster of diameter 
larger than c log n would have to cross one of these parallelograms "in the easy 
direction", which - using the exponential decay property of percolation - has a 
probability at most 

C t 2 e- Cip " )closn < C n 6 - Cip " )c , (3.3) 
that tends to as n — * oo for c large enough. 
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In fact, the case t > n 3 can be handled more directly. We dominate our 
model by using N ~ V(2n): N > n with high probability, and we have 

n(z) = n t (\\z\\)[l + 0(1A 3/4 )] = 0(l/n 3 ), 

so that P2n,t{z) = 0(1 /n 2 ). Hence, the probability to observe two neighboring 
sites is at most 

V" ^ P2n,t( z )P2n,t( z ') < 6 max P2n,t( z ') X ( P2n,t( z ) I = °( ~ )) 

— ' — z'EV T \ — / \ Ti / 

z€V T z'~z \zev T / v 7 

since ^ zE vt P2n,t( z ) ^ 2n - 

Case (HI]): e (log " )Q < t n < X'n for some a E (0, 1) and A' < A c 

Let us turn now to case (p}. First of all, the hypothesis on t n ensures that 

CyVi < r* it < C a (logt)^ 2a Vt (3.4) 

Let us first study the behavior of p n .t around r* t : we have 



-Q^Pn,t{r) = (l-p„,t(r)) x (log(l -p„,t(r))) x f — J. 



(3.5) 



The hypothesis t n < X'n also implies that 

Pn,t (0) > 1/2 + 2(5 (3.6) 

for some <5 > (like in Eq. (|3.ip ). We can thus define 

r- >t = {p n>t )- l {l/2 + 5) and r+ t = (p n ,t)-\l/2 - 5) (3.7) 

as we defined r* t) so that p n ,t(r) is bounded away from and 1 for r S [r~ t , r„ t ]. 
Note that it also implies that nTt t (r) is bounded away from and oo. A formula 
similar to Eq. (|2.16p apply in this case, in particular we have 

C±V~t < < t < C^logty/^Vt. (3.8) 

Hence, we get from Eq. (|3.5|l that for r £ \Tnu r nt\i 

Ci(log*) 1/2a . d C 2 
~ d-r PnAr) ~ - VT (3 ' 9) 

We deduce 

2 ^ ( r - r n, f ) <i?n,t(r) < - - -^(r-r„ it ) (3.10) 
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for r e [r* ( ,r„ t ], and the same is true with inequalities reversed for r e 
[r~ t ,r* t ]. Now, for any z in the annulus S r - r + , 

PnM ~ PnAz) = e-" ff * (l|2|l) (l ~ e»(**(ll*ll)-**W)), (3.11) 

and using Eq.jH]), we get n(7f t (||z||) - Tr t (z)) = (mr 4 (||z||)) x 0(l/i 3 / 4 ) = 
0(l/t 3 / 4 ), so that 

Pn,t{z) ~ Pn,t{z) = 0(1A 3 / 4 ). (3.12) 

Hence, Eq . (|3. 10|1 is also valid with p n< t(z), taking different constants if necessary: 

PnM = 1-M(|| 2f ||- < >t ), (3.13) 

where C( < e(z) < C^logi) 1 / 2 ". 

We now apply a Poissonian approximation, as explained in Section [231 m 
a well-chosen subset A of sites - typically an annulus of width t 2 / 1+t around 
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dS r * j but we need to enhance it a little bit. More precisely, we consider 
the configuration depicted on Figure [2 we take as a subset A of sites the 
annulus S r * t -2t 2 / 7 +>,r* t +2t 2 / 7 +«j together with the smaller annulus and 
some strips connecting these annuli, s^f = \r~ t /2, r* t -0.99 x i 2 / 7+c ] x [0,t 2 / 7+£ ], 
S3 = [0,r~ t ] x [0,i e ] and some vertical strip in the middle. We also add 
external strips, s+ = [r* t + 0.99 x t 2 / 7+£ , 2r+ t ] x [0, t 2 / 7+t ] and s+ = [r+ t ,t] x 
[0, t e ], connected by a vertical strip sj. 

We claim that for this choice of A, we have 

TTt(A) = 0(t- 3/u+2e ). (3.14) 

Indeed, this is a direct consequence of 7r t (z) = 0(l/t) and r* n t — 0(^tlogt), 
except for the most external strip = [r£ t ,t] x [0,i £ ]: for this strip, we have 
(recall that r+ t > CjV*) 

( t 9/16 V 
f J2 \ e ^ jt +0(i e xt xe- tl/S / 2 ) 
r =C+ x/t ' 
= 0(t- 1 /2+e )i 

using Eqs.pTS]) and ([273]) . 

In both annuli 5 r * t _2t 2 / 7 + e .r* t -t 2 / 7 +« and S r * t+t 2/7+ e j ,» t +2t 2 / 7 +«i the corre- 
lation length is at most 

L (l ± ^t 2/7+£ ) = ± C 2 <- 3 / 14+e ^) « t 2 /W3 ; (3 . 15) 

the property of exponential decay with respect to L (Eq. (|2.9|l ) thus implies that 
there are circuits in these annuli, an occupied circuit in the internal one and 
a vacant circuit in the external one: indeed, these circuits can be constructed 
with roughly r* t /t 2 ^ 7 sa i 3 / 14 overlapping parallelograms of size 2t 2 / 7 x t 2 ! 7 
(this size is at the same time » ^2/7-4e/3^ g0 tnat crossm g probabilities tend 
to 1 sub-exponentially fast, and <C t 2 / 7+e ). For the same reason, there are also 
occupied crossings in and s^~, and vacant crossings in and sf - Since the 
parameter p n .t is at least 1/2 + S in S r - , there is also an occupied crossing 

in S3 , and an occupied circuit in St^2t" ■ Since p n> t is at most 1/2 — 8 outside 
of S r + , there is also a vacant crossing in S3 . As a consequence, the potential 
interfaces have to be located in the annulus S r - t -2t 2 / 7 +> ,r* t +2t 2 / 7 +>- 

Let us now consider of order r* t /t 2 ' 7+e ~ ^/u-e su b-sections of length 
t 2 / 7+e in the main annulus as on Figure 03 If we use the same construction 
as for uniqueness in Proposition 7 of [15] (see also p~7]), we can then show 
that with probability tending to 1, there is an edge e (on the dual lattice) in 
S r » _ { 2/7+« ir » t+t 2/7+e with two arms, one occupied going to dS r , t _ 2 t 2 / 7 +> and 
one vacant to dS r » +2t 2/7+«. In this case, the interface is unique, and can be 
characterized as the set of edges possessing two such arms. 
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Figure 3: By considering sub-sections of the main annulus, we can show the 
existence of an edge in S r * _ t 2/7+ c>r * t+t 2/7+ e possessing two arms as depicted, 
which implies uniqueness of the interface. 

To see that the fluctuations of the interface are larger than t 2 ^ 7 ~ e , we can 
use "blocking vertical crossings" as for Theorem 6 in [15]: the definition of L 
indeed implies that in each rhombus of size 2< 2 / 7_e centered on dS r * , there is 
a "vertical" crossing with probability bounded below by 1/4. 

We have just seen that the interface can be characterized as the set of two- 
arm edges. This allows to describe its fine local structure, by comparing it to 
critical percolation interfaces like in [15]. Hence, its length - number of edges - 
behaves in expectation as 

E[L n ] « e^Vtit 2 ' 7 )- 1 '* = t 5 ' 7 , (3.16) 

«2 = 1/4 being the so-called "two-arm exponent" measuring the probability 
for critical percolation to observe two arms: the probability to observe two 
such arms going to distance n decreases as n~" 2 . Moreover, a decorrelation 
property between the different sub-sections also holds, as in [15) . so that L n is 
concentrated around its expectation: more precisely, L n /E[L n ] — > 1 in L 2 as 
n — ► oo. 

Remark 2. We could also construct scaling limits as in f!7\j . using the tech- 
nology of Aizenman and Burchard fTj to prove tighness. One can for example 
look at the portion of the interface close to the x-axis and scale it by a quantity 
of order t 2 / 7 (the right quantity is called "characteristic length for the gradient 
percolation model" in '17]). The curves we get in this way are of Hausdorff 
dimension 7/4, as for critical percolation. 

Remark 3. Here case contains simultaneously the two cases t ~ An (for 
A < \ c ) and t ~ n a . However, note that these two regimes can be distinguished 
by the size of the 'transition window", where the parameter decreases from 1/2+5 
to 1/2 — 5, as shown on Figure^ 

• When t ~ An, the transition takes place gradually with respect to \fi. The 
parameter at the origin is bounded away from 1, so that there are "holes" 
everywhere, even near the origin. 

• When t ~ n Q , the transition is concentrated in a window of size ~ \/i/y/logt 
around r* X y/TTogi. The parameter at the origin tends to 1, and the 
macroscopic cluster is very dense around the origin. 
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Here we used that around r* t , the local behavior of p n j is similar for all t such 
that e^ 08 ™* 1 < t < X'n, with at most an extra (logf) 1 / 2 " factor that does not 
affect the exponents. 



4 Model with a source 

We briefly discuss in this section a dynamical model with a source of particles 
at the origin (think for instance of an ink stain spreading out). We start with 
no particles, and at each time t > 0, new particles are added at the origin with 
rate \i: the number n t of new particles at time t has distribution for some 
fixed [i > 0. Once arrived, the particles perform independent simple random 
walks. As shown on Figure OH we observe in this setting a macroscopic cluster 
around the origin, whose size increases as the number of particles gets larger 
and larger. 

The number of particles N t {z) on site z at time t is given by a sum of 
independent Poisson random variables: 

N t (z) ~ V(jin t (z)) + ...+ V(fin (z)) ~ V( m (z)), (4.1) 

with 

t 

p t (z) = J2'*u{z)- ( 4 - 2 ) 

u=0 

We thus introduce as for the previous model 

q li ,t(z) = l-e-^'\ (4.3) 

For any value of fx, this parameter tends to 1 near the origin, and this 
dynamical model turns out to behave like the previous model in dense phase. 
We show that the constructions of the previous section can easily be adapted 
to this new setting. 

Lemma 4. For each e > 0, there exist constants C'i = Ci{e) such that 

p t {z)>C 1 \ogt-C 2 (4.4) 
uniformly on sites z such that ||z|| < i 1 / 2-6 . 

Hence in particular, q^j{z) — > 1 uniformly on S t i/2-e. 
Proof. We write 



Since \\z\\ 2 /u < t^/u < 1, and E*^-,. 1/u 7 / 4 < J2u=i l l uVi < °°> we S et 
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(a) t = 10. 



(b) t = 100. 



(c) t = 1000. 



Figure 5: A "stain" growing as the time t increases, with new particles arriving 
at rate [i = 50. 



whence the result, since 
t 

a 



7, = ( lo §* + 7 + °W) - ((! - 2e ) iog^ + 7 + o(l)) 

M = t 1 - 2e 

= 2elogi + o(l). 



Let us now consider for r > 



Lemma 5. We have 



uniformly on sites z such that \\z\\ >t 7 / 16 . 
Proof. We write 

t 3 - 74 t 

M=0 U=t 3 / 4 + l 

,3/4 



M=0 U=*3/4_|_l \ V 



□ 



Pt(r) = ^ / — (4.6) 



2tt J r 2 1 1 u 
and also 

^, t (r) = l-e-^W. (4.7) 



/>'(-)=/"( II II ) + ^(^97i«) ( 4 -8) 
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For each < u < i 3 / 4 , we have by Eq. lpTg]) 

n u (z) < Ce-^l^ < Ce- tl/8 / 2 



(4.9) 



since 1 1 z \ \ 2 > i 7 / 8 , so that the corresponding sum is a 
We also have 

i f / 1 — 

U 7 / 4 



«=*3/4 + l 

Let us now estimate the main sum. We can write 
t , t 
-e 



t 9/m 



u=t 3 / 4 + l u =i 3 / 4 + l 

with f\(x) = er x l' x jx and A = ||z|| 2 /i > i -1 / 8 , and we have 

t nu/t 



(4.10) 



(x)dx 



=t 3 / 4 + l 

By an easy calculation, 



< J2 



/a( y J -/a(x) 



dx. 
(4.11) 



^) = i("-l)e-^=^(z 



J / A 

a; 



with g(y) = y 2 (y — l)e a , that is bounded on [0, +oo). Hence for some M > 0, 



so that for x <E [(u — l)/£, w/t], 



/a( - )-fx(x) 



< c't 1 / 4 



u 

X 

t 



< 



a 

43/4' 



(4.12) 



and finally 



E 

U=t 3 / 4 + l 



u/t 
(«-!)/< 



Putting everything together, we get 



Making the change of variable u — \\z\\ 2 /tx, we obtain 
Pt{z) = 



2 A 



/ t3/4 e-« / 1 
du + O „ , 1K 



(4.13) 



(4.14) 
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which allows us to conclude, since 

f °° — du = o( -^L-NI 2 /* 3 ^ = OCt-V'e"^). (4.15) 

□ 

Since ?^,t(r) is strictly decreasing, tends to as r — > oo, and to 1 as r — > + , 
we can introduce 

< t = (f^) -1 ^), (4.16) 
and an easy calculation shows in this case that 




(4.17) 



with F{x) = J °° ^——du, that is a decreasing bijection from (0, +oo) onto itself. 
We also consider 

r^t = (%,t)- 1 m) and r+ t = (^^(l/*). (4.18) 

and we have 

r% = C^Vt. (4.19) 
If we compute the derivative of g M ,t(r) with respect to r, 

9 _ _ , V3 M e- r2 / t 2r 

^.*« = (l-^.«(r))x — X-j^Xp 

we get, similarly to Eq. (|3.10p . that 

yr)^-^( r -r; f ), (4.20) 

for r e [r~j,rt ( ], where Ci < e(r) < C2. We are thus in a position to apply 
the same construction as in the previous section (Figure [2]), and the properties 
of TheoremQ] (case (p}) still hold: for any small e > 0, with probability tending 

to 1 as t — » 00, 

(a) there is a unique interface 74 surrounding dSf , 

(b) this interface remains localized in the annulus S r * _ t 2/7+e )T .* t+t 2/7+ 5 around 
dS r * t , and its fluctuations compared to the circle r — r* t are larger than 
i 2 / 7-e , both inward and outward, 

(c) its discrete length L t satisfies i 5 / 7_e < L t < t 5 ^ r+e . 
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Remark 6. We could also modify our model by adding a fixed number hq of 
particles at each time t: the same results apply in this case, by using a Poisso- 
nian approximation in the same set A of sites as in the previous section. Indeed, 
we can ensure that the configuration in A coincides with the "Poissonian con- 
figuration" with probability at least 1 — pt{A), for 
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